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A flux liquid can condense into a smectic crystal in a pure layered superconductor with the magnetic 
field oriented nearly parallel to the layers. Similar order can arise in low temperature 4 He films 
with a highly anisotropic periodic substrate potential. If the smectic order is commensurate with 
the layering, this periodic array is stable to quenched random point-like disorder. By tilting and 
adjusting the magnitude of the applied field, both incommensurate and tilted smectic and crystalline 
phases are found for vortex arrays. Related variations are possible by changing the chemical potential 
in the helium system. We discuss transport near the second order smectic freezing transition, 
and show that permeation modes in superconductors lead to a small non-zero resistivity and a 
large but finite tilt modulus in the smectic crystal. In helium films, the theory predicts a nonzero 
superfluid density and propagating third sound modes, showing that the quantum smectic is always 
simultaneously crystalline and superfluid. 

PACS: 74.60.Ge,74.40.+k 



I. INTRODUCTION 

The discovery of high temperature superconductors, 
with their broad fluctuation regime, has emphasized the 
inadequacy of the. conventional mean-field description of 
critical behavior .El Early attempts to apply field theoretic 
methods to the Ginzburg-Landau (GL) free energy, iow- 
ever, are of limited applicability in low dimensions. EfB 

Instead, superconducting fluctuations in low dimen- 
sions are now understood in terms of vortices, which 
emerge as the low energy degrees of freedom of the 
Ginzburg-Landau theory. The phases of the supercon- 
ductor within this picture are analogous to states of con- 
ventional matter, except that they are composed of flux 
lines instead of molecules. In fact, because the vortices 
are extended objects, the system most closely resembles 
a collection of quantum bosons .□ As the thickness of the 
superconductor approaches infinity, the effective "tem- 
perature" of this bosonic system goes to zero, and inter- 
esting strongly-correlated phases can emerge. 

One difference between the flux line array and true as- 
semblies of bosons is that in the former, the effects of 
the embedding medium are more often dramatic. In- 
deed, without careful preparation, most high tempera- 
ture superconducting samples are dominated by internal 
defects which tend to disorder the vortex array. Only at 
reasonably high temperatures, when the fluxons are best 
described as a liquid, can the effects of these random pin- 
ning centers be neglected.El At lower temperatures, disor- 
der may induce subtle types of glassy order ,ETQ or simply 
force the system to remain a liquid with very sluggish 
dynamics .Ej 

The layered structure of the copper-oxide materials it- 
self provides a non-random source of pinning S3 At low 



temperatures, the c-axis coherence length £ c o ~ AA < 
s ~ 12A, where s is the lattice constant in this direc- 
tion. Vortex lines oriented in the ab plane are attracted 
to the regions of low condensate electron density between 
the Cu02 layers. Such a periodic potential for true two- 
dimensional bosons could be induced by an anisotropi- 
cally corrugated substrate, possibly leading to the obser- 
vation of the effects described here in He 4 films. 

Previous work on intrinsically pinned vortices has fo- 
cused on the low temperature fluctuationless regime, in 
which the vortices form a pinned elastic solid. Near T c , 
however, when thermal fluctuations are important, en- 
tirely different phases can exist. These thermally fluc- 
tuating states are particularly interesting experimentally 
because hystcrctic effects arc weak and equilibrium trans- 
port measurements are more easily performed than at low 
temperatures. Our research is motivated by the recent 
experimental work of Kwok et. al.E3 who observed a con- 
tinuous resistive transition in YBa2Cu307 for fields very 
closely aligned (0 < 1°) to the ab plane. A preliminary 
version of our results appeared in Ref . 

To explain the experiments, the interplay between 
inter-vortex interactions and thermal fluctuations must 
be taken into account in an essential way. The experi- 
ments of Ref. |lj seem to rule out conventional freezing, 
which is first order in all known three-dimensional cases. 
The additional observation of a strong first order freez- 
ing transition for 8 > 1° suggests that point disorder 
is relatively unimportant at these elevated temperatures 
(strong point disorder would destroy a first order freezing 
transition). In addition, an attempted fit of the data to 
a dynamical scaling form yieldecLexponents inconsistent 
with vortex or Bose glass valuesJlj Instead, we postulate 
freezing into an intermediate "smectic" phase between 
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the high temperature flux liquid and a low temperature 
crystal/glass. Such smectic freezing, as discussed by de 
Gennes for the nematic-smectic A transition,E-3 can oc- 
cur via a continuous transition in three dimensions. The 
vortex smectic state is richer than its liquid crystal coun- 
terpart, however, for two reasons. First, the existence of 
a periodic embedding medium (i.e. the crystal lattice) in 
the former leads to eammensurability effects not present 
in the liquid crystal.Ej In addition, the connectedness of 
flux lines leads to constraints with no analog for point- 
like molecules. As we show below, the onset of smectic 
order should be accompanied by a steep drop in the re- 
sistivity and a rapid increase in the tilt modulus for fields 
which attemp to tip vortices out of the Cu02 planes - 
see Fig.|Io|. 

FIG. 1. Schematic structure functions in the (a) liquid, (b) 
smectic, and (c) solid phases. For simplicity, we have illus- 
trated the case for a square lattice. 

The smectic phase may also be distinguished experi- 
mentally using neutron scattering, which measures the 
Fourier transform of the magnetic field two-point corre- 
lation function (see Fig|lJ). We assume a magnetic field 
along the y-axis and Cu02 layers perpendicular to z. 
The vortex liquid structure function shows the usual dif- 
fuse liquid rings, as well as delta-function Bragg peaks 
at q z — 27m/ s (for integral n), representing the "im- 
posed" vortex density oscillations from the Cu02 layers. 
On passing to the smectic state, additional peaks develop 
at wavevectors q z = 27m/ a, interlacing between those al- 
ready present in the liquid. The new peaks represent 
the broken symmetry associated with preferential occu- 
pation of a periodic subset of the layers occupied by the 
vortices in the liquid. At lower temperatures in the vor- 
tex solid, further peaks form for q x ^ 0, producing the 
full reciprocal lattice of a two-dimensional crystal. 

Our analysis leads to the phase diagrams shown in 
Fig.^. Upon lowering the temperature for H c = and a 
commensurate value of Hb, the vortex liquid (L) freezes 
first at T s into the pinned smectic (S) state, followed by 
a second freezing transition at lower temperatures into 
the true vortex crystal (X). When H c ^ 0, tilted smec- 
tic (TS) and crystal (TX) phases appear. The TS-L and 
TX-TL transitions are XY-like, while the TS-S and TX- 
X phase boundaries-are commensurate-incommensurate 
transitions (CITs)£3 At larger tilts, the TX-TS and TS- 
L phase boundaries merge into a single first order melting 
line. As ifj, is changed, incommensurate smectic (IS) and 
crystal (IX) phases appear, again separated by CITs from 
the pinned phases, and an XY transition between the IS 
and L states. If the low— Hb commensurate smectic (S) 
phase corresponds to, say, 5 Cu02 plane periodicities per 
vortex layer, the high— iff, S phase will represent a state 
with 4 Cu02 plane periodicities per sheet. We show that 
the commensurate smectic order along the c axis is stable 
to weak point disorder, in striking contrast to the trian- 
gular flux -Lattice which appears for fields aligned with 
the c axis£3 . This stability should increase the range 
of smectic behavior relative to the (unstable) crystalline 
phases when strong point disorder is present. 



FIG. 2. Phase diagrams in the (a) H c -T, and (b) H b -T 
planes for the disorder-free vortex system. The fuzzy lines 
indicate first order transitions. In (b) two possible topologies 
are shown connecting different commensurate states. 

The high temperature flux liquid (for fields in the ab 
plane) is of some interest in its own right. Consider first 
one vortex line wandering along the y-axis, as shown 
schematically in FigJ^. This line is subject only to ther- 
mal fluctuations and a periodic pinning potential along 
the z-axis, provided by the Cu02 planesO If thermal 
fluctuations are ignored, the vortex acts like a rigid rod_ 
and will be localized in one of the potential minima.E£l 
This localization assumption, however, is always incor- 
rect in the presence of thermal fluctuations, provided 
the sample is sufficiently large in the y-direction. As 
L y — > oo, the statistical mechanics of this single wander- 
ing line random walking in directions perpendicular to y 
leads inevitably to equal probabilities that the vortex is 
in any of the many possible minima along z. 

On a more formal level, this probability distribution 
P(z) is given by the square of the ground state wave 
function of the Schrodinger equation in a periodic poten- 
tial - see section ID below. The jumps shown in Fig.|| 
across Cu02 planes are represented by quantum mechan- 
ical tunneling in imaginary time. According to Bloch's 
theorem, this tunneling leads to P(z) = \ipk=o{z)\ 2 , 
where the Bloch states in general have the form ipk(z) — 
exp(ikz)u(z), with u{z) a function with the periodicity 
of the pinning potential. The resulting probability distri- 
bution is shown schematically on the right side of Fig|| 

FIG. 3. Wandering of a single flux line (solid curve), lead- 
ing to an extended probability distribution P(z) given by a 
k = Bloch wavefunction. Other vortex trajectories (repre- 
sented by the dashed curve) will generate similar probability 
distributions, unless interactions lead to crystalline or smectic 
order. 

Now suppose an additional line is added to the sys- 
tem. As suggested by the trajectory of the dashed curve 
in Fig.[|, it too will wander from plane to plane. Although 
the two flux lines interact repulsively, they can wander 
and still avoid each other by using the ^-coordinate or by 
never occupying the same minimum at the same value of 
the "imaginary time" coordinate y. Thus both flux lines 
generate a delocalized probability distribution and oc- 
cupy the same k = Bloch state. At high temperatures 
or when the lines are dilute, we expect for similar rea- 
sons macroscopic occupation of the k = Bloch state 
in the equivalent boson many body quantum mechan- 
ics problem, similar to Bose-Einstein condensation. In 
this sense, the flux liquid is indeed a "superfluid." The 
presence of numerous "kinks" in the vortex trajectories 
insures a large tilt response for fields along z and a large 
resistivity for currents along x. The various symmetry- 
breaking crystalline or smectic states which appear at 
low temperatures or higher densities arise because of the 
localizing tendency of the interactions. The density of 
kinks is greatly reduced in these phases. 

The remainder of the paper is organized as follows. In 
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section |l[ several models are introduced which will be 
used to analyze the layered superconductor. Sections III 
and IV discuss the effect of intrinsic pinning on the liq- 
uid and crystal phases, respectively, and show how smec- 
tic ordering is encouraged on approaching the interme- 
diate regime from these two limits. A Landau theory 
for the liquid-smectic transition is introduced in section 
0, and the critical behavior is determined within this 
model. The nature of the commensurate smectic phase 
itself is explored in section VI, through a computation 
of the response functions. In section VII, it is shown to 
have "supersolid" order, similar to the supersolid crystal 
phase recently proposed at high magnetic fields along the 
c axis.ta The additional phases which arise f or lar ge in- 
commensurate fields are described in section VIII. Sec- 



tion [X details the modifications of the phase diagram 
when weak point disorder is present, and, in particular, 
demonstrates the stability of the smectic state. Conclud- 
ing remarks and the implications of these results for he- 
lium films on periodically ruled substrates are presented 
in section [x|. 



II. MODELS 

At a fundamental level (within condensed matter 
physics), layered superconductors may be modeled as a 
positively charged ionic background and a collection of 
conduction electrons, which can pair via the exchange 
of phonons, excitons, magnons, etc. Since such a mi- 
croscopic theory of high temperature superconductors is 
lacking, we must resort to more phenomenological meth- 
ods. There are, nevertheless, a variety of differing levels 
of description available, several of which will be used in 
the remainder of this paper. 



A. Static Models 

The most basic of our models is the familiar Ginzburg- 
Landau theory, which is an expansion of the free energy 
of the superconductor in powers of the order parameter 
field * GL , 
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where D M = 9 M — 2%A^/(po, <fio = hc/2e is the flux quan- 
tum. Here H is the applied external magnetic field. To 
establish notation, we choose x, y, and z respectively 
along the a, b, and c axes of the underlying cuprate crys- 
tal. The diagonal components of the effective mass ten- 
sor are then m x = m y = m and m z — M . For conve- 
nience in what follows, we also define the anisotropy ratio 

7 = y/m/M = Xab/Xc = £,c/£,ab < 1. 



Eq.|l| provides a powerful means of understand- 
ing superconducting behavior, including the effects of 
anisotropy. In layered materials, however, the theory 
must be modified to allow for coupling of the supercon- 
ducting order to the crystalline lattice. One such de- 
scription, in which the superconductor is regarded as 
a stack of Josephson-coupled layers, is the Lawrence- 
Doniach model. For our purposes, however, it is sufficient 
to consider a "soft" model for the lattice effects, in which 
the coupling a is allowed to be a periodic function of z 
with period s equal to the periodicity of the copper-oxide 
planes. 

Because the high temperature superconductors are 
strongly type II, it is appropriate to use London theory 
over a large range of the phase diagram. In this limit, 
variations of the magnitude of the order parameter are 
confined to a narrow region within the core of each vor- 
tex. Because the resulting London equations are linear, 
a complete solution can be obtained— for the free energy 
of an arbitrary vortex configuration.Ell For our purposes, 
it is sufficient to consider an approximate form in which 
the tilt moduli are local-and the interactions between 
vortices occur at equal 
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dyV[v ±i (y)-r ±j (y)], (2) 



where V[r x ] » 2 1 e K Q ( A /x 2 + z 2 /j 2 /\ ab ) {K a {r) is a 
modified Bessel function), eo = (<f>o/ '47rA a f,) 2 ), and the 
stiffness constants obtained from anisotropic Ginzburg- 
Landau (GL) theory are eii = eo7 and 6^ = 60/7 (see, e.g. 
Ref. Vp[z] is a periodic potential (Vp[z + s] = Vp[,z]) 
taking into account the effects of the layering. 

Either Eq.[j] or Eq.|| may be used at finite temperature, 
by calculating the partition function 



Z = Tre- F / fcBT , 



(3) 



where the trace is a functional integral over \&gl or 
the set of vortex trajectories {1^(2)}, for the Ginzburg- 
Landau and London limits, respectively. 

In the London case, this trace may be formally per- 
formed by recognizing Eq.[| as mathematically identical 
to the Feynmann path integral for the first quantized 
imaginary time Green's function of interacting bosons. 
Within this boson analogyu the Green's function may 
also be calculated using a coherent state path integral 
representation. The "action" for these bosons is 



S hoson = f d 3 r^(Td y - ^d 2 z - ^d 2 x 
J Zej_ zeii 



d 3 rV P (z)n(r) 



+ d 2 r ± d 2 r' ± dy-V(r ± -r' ± )n(r ± ,y)n(r' ± ,y). (4) 



where if} is the complex coherent state boson field, and 
n(r) = V^( r )V'( r )- It is convenient to rescale x — > 



3 



(e^/e^Y^x and z — > (e r ||/e^) 1 / 2 z to obtain the isotropic 
Laplacian = d 2 + d 2 . Eq.[| becomes 



^(Td y -^Vi-^-V P (z)n(r) 



+ ^ j d 2 r ± d 2 r' x dyV(r x -r' x )n(r x ,y)n(r' x ,y), (5) 



with e = */ej|eI and the rescaled potentials Vp(z) = 
Vp(z7 2 ) and V(r±) = V(x/j 2 ,z-f 2 ). The action is used 
via 



[# f ][d# 
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(G) 



to calculate the grand canonical partition function Z at 
chemical potential /x per unit length of vortex line. EqsJ^ 
and p| may also be obtained dixectly froma limiting case 
of Eq|l] via a duality mapping Bj 



B. Dynamical Models 

Calculation of dynamical response functions such as 
the resistivity requires a model for the time dependence 
of the superconductor. We will do this within the Lon- 
don framework, treating the vortex lines as the dynam- 
ical degrees of freedom. In the overdamped limit, the 
appropriate equation of motion is then 



rrxi(2/) 
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where T is a damping constant, and f; is the force on the 
i th flux line, including both external forces and random 
thermal noise. 

Eq.fj] is useful in describing the properties of small 
numbers of vortices. To understand the bulk behavior 
of dense phases, however, we need an extensive theory. 
Such a model for the liquid phase was constructed on 
physical and symmetry grounds in Ref . ^4] in the hydro- 
dynamic limit. Because we intend to go beyond simple 
linearized hydrodynamics, we require some knowledge of 
the non-linear form of the bulk equations of motion. 

We proceed by first defining the hydrodynamic fields 



n(r) =^<5[r_L -rxi(y)], 

i 

i d y 



(8) 

(9) 
(10) 



Conservation of magnetic flux is embodied in the con- 
straint 

d v n + V± • t = 0. (11) 
Eq.|| may be rewritten as a function of n and r, formally 

^London [Wi/)}] = F[n(r) , T (r)}. (12) 



Specific forms for F[n, t] will be used as needed. 

Eq.Q completely specifies the dynamics of n and t. 
Differentiating Eqs.|| and || then leads immediately to 
hydrodynamic equations of motion. Details are given in 
appendix [A|. One finds 

d t n + Vi • j„ = 0, (13) 

d t T a + dpjpa = dyjv,a, (14) 

where the density and tangent currents are 
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(17) 



Because the constitutive relation for the tangent cur- 
rent includes t' 37 , Eqs.|l3|-{i"7j do not form a closed set. 
Vortex hydrodynamics, however, leads us to expect that 
n and r provide a complete long-wavelength description 
of the system. We therefore adopt the truncation scheme 
(i.e. t^ 1 is preaveraged at equilibrium). 
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Averaging via Eq.^| gives 



^ -> k B Tn 



ex 



(18) 



To complete the dynamical description, we must specify 
r and f . Matching to the liquid hydrodynamics of Ref. 
|24| relates V to the Bardeen-Stephen friction coefficient 
7bs = ^or, and gives the driving force 



f =^JAy + J7, 

c 



(19) 



Here J is the applied transport current density and r\ (r) 
is a random thermal noise. 



III. INTRINSIC PINNING IN THE VORTEX 
LIQUID 

To better understand the interplay of thermal fluctu- 
ations and layering of the superconductor, it is useful to 
first consider the behavior at high temperatures in the 
liquid state. As discussed by Marchetti,c3 the physics of 
a single vortex line in the liquid is well described by a 
hydrodynamic coupling to the motion of other vortices. 
In the dense limit, this medium is approximately uni- 
form, and does not significantly effect the wandering of 
an isolated vortex. To estimate the effects of intrinsic 
pinning, it is thus appropriate to consider a single vortex 
line, oriented along the a — b plane. 

For T > 80K, £ c « £ c0 (l - T/T c )-^ 2 > s. In this 
limit, the copper-oxide planes act as a smooth periodic 
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potential on the vortex. The magnitude of this potential 
per unit length isE£l 



■ > \ 10 2 e 7 ( — 



5/2 



-15.8£ e /a 



For a single vortex, Eq.0 reduces to 



dy{ 



dx 
dy 



V P (z) 



(20) 



(21) 



The periodic potential Vp[z] — U p f p [z/s] 7 where f p (u) = 
f p (u + 1) is a smooth periodic function with magnitude 
of order unity. 

The x-displacement decouples in Eq|2l], and may be 
integrated out to yield 



([x(y) - -t(O)] 2 ) ~ D x y, 



(22) 



with D x = k^,T/e\\. The z-dependent part of Eq|2l] is 
identical to the Euclidean action of a quantum particle 
of mass £j_ in a one-dimensional periodic potential Vp(z), 
with y playing the role of imaginary time. The single flux 
line partition function, 



[dz(y)]e 



-F v [z]/k B T 



(23) 



with fixed endpoints, maps to the Euclidean Green's 
function for the particle, with k B T replacing Ti. 

In the quantum-mechanical analogy, the particle tun- 
nels between adjacent minima of the pinning potential, 
leading, as discussed in the Introduction, completely de- 
localized Bloch wavefunctions even for extremely strong 
pinning. The "time" required for this tunneling maps 
to the distance, £kink, in the y-direction between kinks 
in which the vortex jumps across one Cu02 layer. The 
WKB approximation gives 



L 



kink 




Ups/knT 



(24) 



Eq.^3] may also be obtained from simple scaling con- 
siderations. The energy of an optimal kink is found by 
minimizing 



h 



w 



UpW 



(25) 



over the width (in the y direction) w, giving w* ~ 
\feJjTTpS and /* ~ y/exU p s. Such a kink occurs with a 
probability proportional to exp(— fl/k^T) in a length w. 
The condition i kink /wexp(— fi/k B T) ~ 0(1) then gives 

Eq.|§ 

When the sample is larger than Lkink along the field 
axis, the flux line will wander as a function of y, with 

([z(y) - z(0)} 2 ) ~D z y, (26) 

where the "diffusion constant" D z w s 2 /£kink- 

For y/e±U p s < k B T, the pinning is extremely weak, 
and the WKB approximation is no longer valid. Instead, 



the diffusion constant D z « k B T/e±, as obtained from 
Eq.^T] with Up — 0. At much lower temperatures, when 
£ c <C s, the energy in Eq.|2(] must-be replaced by the 
cost of creating a "pancake" vortexED between the Cu02 
planes. In this regime, L kink ~ ^ab(s/£, c ) e ° s/kBT ■ 

gr T w 90K, as in the experiments of Kwok et. 
S, c /s w 2.3, and Eq.|^ gives y / e±U p s/k B T < 1, 
indicative of weakly pinned vortices in the liquid state. 
The transverse wandering in this anisotropic liquid is de- 
scribed by a boson "wavefunction" with support over an 
elliptical region of area k B TL y l ^J^t± with aspect ratio 
Ax/Az = 7 _1 w 5 for YBa2Cu3C>7. For L y i=s 1mm, a 
typical sample dimension along y, the dimensions of this 
ellipse are of order microns. Since typical vortex spacings 
at the fields used in Ref. |l2| are of order 400A, these flux 
lines are highly entangled. 

To understand the bulk propertices of the vortex liq- 
uid, it is usefu l to employ the hydrodynamic description 
of section [II Bj Jr^ the liquid, the appropriate form of the 
free energyisOCJ 

FL = hj (^{ c ii(q)IMq)| 2 +c 4 4,||(q)|r,(q)| 2 
+c 44 ,±(q)Mq)| 2 j - Jd 3 rV P [z}Sn(r), (27) 

where Sn = n — no, with uq — B y /<pQ the mean density. 
Here the compression modulus en and tilt moduli 044.^ 
and C44 || are regular functions of q with finite values at 
q = 0. ' 

On physical grounds, we expect intrinsic pinning to 
enter Eq.^ both through an increase in C44 j j_, which 
decreases fluctuations perpendicular to the layers, and 
through the Vp term which tends to localize the vortices 
near the minima in the periodic potential. 

The former effect only acts to increase the anisotropy 
of the liquid. The latter term, however, explicitly breaks 
translational symmetry along the z axis, inducing a mod- 
ulation of the vortex density, 



dze 



'■(6n(r)) 



c\i{q z ,q x = q y = 0) 



V P [q z }. (28) 



This modulation corrects the static structure function, 
5(q) = ((5n(q)5n(-q))/(27r) 3 (5 (3) (q = 0), according to 



S(q) = 5 (q) 



[cii(g z )]' 



:\V P [q z }\ 2 (27T) 2 S(q x )S(q y ), (29) 



where So(q) is the static structure function for Vp = 0. 
Since Vp[z] is a periodic function, the correction term 
shows peaks at the discrete reciprocal lattice vectors for 
which q z is an integral multiple of 2tt/sz. 

The situation is somewhat analogous to applying a 
weak uniform field to a paramagnet, inducing a propor- 
tionate magnetization. Unlike the magnetic case, how- 
ever, the layering perturbation leaves a residual transla- 
tional symmetry under shifts z — > z + s. It is the break- 
ing of this discrete group which we will identify with the 
freezing of the vortex liquid. 



5 



IV. THE CRYSTAL PHASE 

Considerable wodc already exists on intrinsic pinning 
in vortex crystals H We review the essential ideas here, 
and discuss its implications for thermal fluctuations at 
low temperatures. 



A. Zero temperature properties 

To study the effect of layering upon the vortex state, 
we first consider the limit of a weak periodic modulation 
of the order parameter along the z-axis. In this case, the 
resulting (zero temperature) configuration is only slightly 
perturbed from the ideal lattice predicted by Ginzburg- 
Landau (or London) theory. The free energy in this case 
may be written in terms of the phonon coordinates u(r), 
as 



F, 



clastic 



d 3 q K, 



a(3ij 



(2tt) 3 



q a q/3Ui (q)uj (-q) + Fjp , (30) 



where i and j = x, y, a and (3 = x, y, z, and the contri- 
bution to the free energy of the layering potential is 



Here K — \ a b/£,ab is the usual Ginzburg-Landau param- 
eter, and /?a ~ 1.16. 

From Eqs.[30| and |33|, it is clear that for these commen- 
surate fields, the ground state is unchanged, i.e. u = 0. 
Away from these fields, however, the fate of the lattice 
is less obvious. Ivlev et. al.El have shown that, for a 
small deviation from a commensurate field, it is energet- 
ically favorable for the vortex lattice to shear in order to 
remain commensurate with the copper-oxide plane spac- 
ing. Because such a distortion requires some additional 
free energy, it will generally be favorable, in addition, for 
the internal magnetic induction B to deviate from the 
applied field H to allow a better fit to the crystal. This 
Meissner-like effect will be discussed in more detail in 
section 
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For strong layering, such as that described by the 
Lawrence-Doniach model, the pinning effects are much 
more pronounced. When £ c <C s, the magnetic field re- 
mains essentially confined between the CuOi layers, and 
the vortex array is thus automatically commensurate at 
all applied fields. Although such strong confinement pfj 
vortices can lead to interesting non-equilibrium states,E2l 
we will confine our discussion to equilibrium. 



■-IP 



dy 



V P [ 



(31) B. Thermal Fluctuations about the commensurate 

state 



In general the elasticity theory is quite complex due to 
the anisotropy and wavevector dependence on the scale 
of A. Rather than work with a specific form of the elastic 
moduli, we will obtain general expressions in terms of an 
unspecified set of iCy(q) = K a pijq a qp. 

To obtain the correct continuum limit of Eq.|3l], we con- 
sider the possible commensurate states of the layers and 
vortex array. The triangular equilibrium lattice in this 
orientation is described by the two lattice vectors ai = 
C 7 ~ 1 z, a 2 = C7"V2(z+\/37 2 x), with C 2 = 20 o /\/3B y . 
Commensurability effects occur when the minimum z- 
displacement between vortices, C/(2j) = (n/m)s, where 
n and m are integers, chosen relatively prime for definitc- 
ness. This gives the commensurate fields 



r>(m,n) 
y 



2V3TS 



2„2 



(32) 



For simplicity we consider here only the integral states 
with m — 1. In this case, Vp[z„ +u] = Vp[u], and we can 
straightforwardly take the continuum limit 



ip 



d 3 vV P [u z (r)]. 



Such an expression may be derived explicitly from the 
Abrikosov solution for fields near -ff C 2,liil in which case 
the pinning potential is 



Vp[u] ~ cos(2ttw z /s), 
Zirc 



(34) 



where 



4 c(H c2 - H) g c7 2 



exp(-8^/s 2 ). (35) 



Thermal fluctuations of the vortex lattice are described 
by the partition function 



Z^ ic = j[du(r)]eM-F^ c /k B T). 



(36) 



In three dimensions, phonon fluctuations are small, 
and expanding the pinning potential around its minimum 
for small u, gives the quadratic free energy 



clastic 



d 3 q ify(q) , , , A 
—-^3 J 2 Ui(q)«j-(-q) + yK(q) 



(37) 



, P L , - 0]. 



where A = — V P [z 

The displacement field fluctuations can be calculated 
from Eq.37 by equipartition, yielding the general result 



dh 



K zz (q) + A 



k B T J BZ (2tt) 3 K xx (q) {K zz (q) + A) - [K xz (q) 



r, (38) 



d\ 



(33) k B T J BZ (2tt)3 ^(q)(^(q) + A) - [^(q)] : 



r, (39) 



where BZ indicates an integral over the Brillouin zone. 

The effect of the periodic potential is thus to uniformly 
decrease the fluctuations of u z at all wavevectors. For 
A > B 2 /X 2 , this decrease is substantial over the entire 
Brillouin zone, and 



(ul) w k B T ( t \ 

JBZ ( 27r ) K xx 



k B T 



uz (27r)»^ M (q) 

d 3 q 1 
bz(2^A- 



(40) 
(41) 
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In stronger fields, for A < B 2 /X 2 , the only the contri- 
butions from q < \/~K/B are strongly suppressed. For 
1 - B/H c2 « 1 
the ratio 




and pq can be combined to give 



within a given layer. In fact, a configuration in which a 
single line hops from layer n to layer n + 1 corresponds to 
a dislocation in layer n paired with an anti-dislocation 
in layer n + 1 , since 



A\ 2 /B 2 



-(H c2 /B-1) 



-8(Ws) 2 



(42) 



At lower fields and temperatures, one expects the mean- 
field estimate above to break down and AX 2 /B 2 to in- 
crease, possibly settling down to a constant value at low 
temperatures. 

For magnetic fields oriented along the c axis, the Lin- 
demann criterion has been-psed to estimate the melting 
point of the vortex latticea by requiring that ( |ui | 2 ) = 



2 „i 



ct.a 



for 



x, z, with a "Lindcmann number" cl 



0.2 — 0.4. As is clear from Eqs.p8|-p9|, once layering is in- 
cluded, the increased stiffness for u z makes the two ratios 
(u 2 }/a 2 and (u 2 }/a 2 unequal. Indeed, the second ratio is 
strongly suppressed relative to the first. Extending the 
Lindcmann criterion to this situation suggests that the 
strains in u x might be alleviated by a partial melting of 
the lattice without affecting the broken symmetry lead- 
ing to the u z displacements. Such a scenario corresponds 
to the unbinding of dislocations with Burger's vectors 
along the x axis. The phase in which these dislocations 
are unbound is the smectic. 



C. Strongly Layered Limit 

To further elucidate the nature of the smectic phase, it 
is helpful to discuss the limit of very strong layering. In 
this case, the vortex lines are almost completely confined 
within the spaces between neighboring CuOi layers. For 
moderate fields, occupied layers will be separated by sep- 
arated by several unoccupied ones, and the interactions 
between vortices in different layers may be considered 
weak. Because of the strong layering, the out-of-plane 
component of the displacement field u z is suppressed, so 
that the free energy of the system may be written to a 
first approximation as 



F 



layers 



d 2 q± 
(2tt)2 



where qj_ = (q x ,q y ) and u n (q_i_) = u x (q^,ns). For a 
qualitative discussion of smectic ordering, it is sufficient 
to take K x and K y independent of q. 

Eq.^3] neglects both inter-layer interactions and hop- 
ping. The former are included perturbatively via the free 
energy 



r 2tt 

F int . = dxdyvi L cos— (u n+1 



(44) 



where i>il is an inter-layer interaction energy and a is the 
lattice spacing in the x direction. The periodic form of 
the interaction is required by the symmetry under lattice 
translations u — > u + a within each layer. 

Once hopping of flux lines between neighboring oc- 
cupied layers is included, u n is no longer single valued 



Vu k ■ dl = a(S k , n - S k ,n+l) 



(45) 



for a contour surrounding the hopping point (see Fig^). 
Such dislocation-antidislocation pairs, which we will re- 
fer to as large kinks, can be created in neighboring layers 
with a dislocation fugacity yd = exp(— E\\^/k B T), where 
the core energy 



En 



y / e±U p ms, 



(46) 



as estimated from Eq.|25| with s — > ms. Note that the 
dislocation and anti-dislocation must have the same x 
and y coordinates, since misalignment is accompanied by 
an energy cost proportional to the extra length of vortex 
between the occupied layers. 

The full theory described by Eqsjl^ and [l5] plus dislo- 
cations can be studied using a perturbative renormaliza- 
tion group (RG) expansion in Ujl and y^, using tech- 
niques developed for the XY model in a symmetry- 
breaking field.E3 For yd = vxt = 0, the Gaussian free en- 
ergy of Eq.^ describes a fixed line of independently fluc- 
tuating vortex layers parameterized by the dimension- 
less ratio y/K x Ky/kBT. To characterize the order along 
this fixed line, we define a translational order parameter 
(characterizing correlations along the x axis) within the 



,th 



layer by summing over vortex lines according to 



p\\(x,y, 



exp(2irix k n) {y)/a), 



k 



(47) 



where x k (y) is the coordinate of the fc th vortex line in 
layer n at a length y along the field direction. The cor- 
relation function Ct,|| {x, y, n) = (pu (x, y, n)p\\ (0, 0, 0)) is 

then evaluated by inserting x k n ^ (y) = ka + u n (ka, y) and 
converting the sum to an integral via J^ fc — > J dx/a. One 
finds 



(K x q 2 x + K z q 2 z )\u n (q ± )\ 2 , (43) q 



T.|| 



[x,y,n) 



K 



K x x 2 + K y y 2 



-Trk B T/Ka 2 



Sn,Q, (48) 



i.e. quasi-long-range order within the planes. 

This fixed line is always unstable either to interlayer 
couplings, to dislocations, or to both perturbations. The 
linear (in y and vjl) RG flows which determine the sta- 
bility are 



dvih 
dl 

dyd 
dl 



= 2- 



2nk B T 

Ko 2 

Ka 2 \ 
2nk B T ) 



VIL, 



yd, 



(49) 
(50) 



where K = ^K x K y: and I — ln(6/a) is the logarithm of 
the coarse-graining length scale b. 

When k B T < Ka 2 /An, dislocations are irrelevant at 
the fixed line, so that y decreases under renormalization. 
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In this regime, however, vji, increases with I, so that in- 
teractions between the layers are important for the large 
distance physics. To study this regime, one may therefore 



expand the cosine of Eq,44 in u n+ i — u n: obtaining a dis- 
crete version the usual three-dimensional elastic theory. 
In this limit, 



Cx ii [x, y, n) ~ Const. 



(51) 



for large \x\, \y\, or \n\. 

At high temperatures, when k^T > Ka 2 /ir, i>il scales 
to zero while the fugacity yd is relevant. Unbound dis- 
locations on long length scales therefore invalidate the 
elastic theory of Eq.^. Following standard arguments, 
the translational correlation function in such an unbound 
vortex plasma becomes exponentially small, i.e. 

Ct,|| (x, y, n) ~ exp(-f/£ x ), (52) 

where r = ^J(K X / K)x 2 + (K y /K)y 2 + xs 2 n 2 with % 
a constant, and £t is a finite translational correlation 
length. 

For temperatures in the intermediate range Ka 2 /4ir < 
k-gT < Ka 2 /ir, both y^ and vn, are relevant operators. 
The eventual nature of the ordering at long distances 
presumably takes one of the two above forms, though the 
critical boundary at which the system loses long-range 
translational order in the x direction is not accessible by 
this method. 

It remains to discuss translational order along the lay- 
ering axis. Such order is characterized by the parameter 

P±(x,y,n) = ^2exp(2iriz { f! l) (y)/a z ), (53) 
k 

where a z is the distance between occupied vortex layers 
(and therefore an integral multiple of the CuC>2 plane 
spacing, i.e. a z — ms). Because we have, by con- 
struction, confined the vortices to these layers, however, 
z^ = na z for every k, and the exponential in Eq.^3] is 
always unity. Both phases described above, regardless of 
the relevance of Uil and y, therefore retain long-range 
order along the z axis. 

Within the strongly layered model, there are still exci- 
tations which can destroy this transverse ordering. These 
are configurations in which a flux line hops out of an oc- 
cupied layer into one of the (m — 1) unoccupied interme- 
diate planes between it and the next occupied layer. Such 
an excursion costs an energy proportional to the length 
of the vortex segment in the unoccupied layer, so that 
only short intermediate segments occur at low temper- 
atures. These out-of-plane hops reduce the amplitude 
(p±), but do not drive it to zero. At very high tempera- 
tures, entropy may counterbalance this energy and drive 
the free energy cost for such intervening vortices nega- 
tive. Once this occurs, translational order will be lost 
along the z axis as well, and the system will be a true 
liquid. Nevertheless, at intermediate temperatures above 
the unbinding transition for dislocations in the u field but 
below the temperature at which infinite vortices enter the 
intermediate copper oxide planes, we expect the system 



to sustain "one-dimensional" long-range order along the 
z axis, i.e. a smectic state. 

FIG. 4. Hopping of a vortex line between neighboring vor- 
tex layers. 



V. CRITICAL BEHAVIOR 

Having established the possibility of a smectic phase 
approaching both from the crystalline and liquid limits, 
we now focus on the critical behavior near the putative 
liquid-smectic transition, using a Landau order param- 
eter theory. A closely related Landau theory which de- 
scribes a low temperature smectic-crystal transition is 
discussed in Appendix [B|. The natural order parameter 
to describe the smectic ordering is p± defined in Eq,53. 
To simplify notation, we define a new field $ = px, so 
that, in the continuum notation (i.e. outside the strongly 
layered limit), 



n(r) « n Q Re {l + $(r)e~ l<?z } , 



(54) 



where no is the background density, and q = 2n/a z is 
the wavevector of the smectic layering. The complex 
translational order parameter ^(r) is assumed to vary 
slowly in space. The superconductor is invariant under 
translations and inversions in x and y, and has a discrete 
translational symmetry under z — > z + s, where s is the 
Cu02 double-layer spacing. From Eq.|54|, these periodic 
translations correspond to the phase shifts <& — * $e _lgs . 
We continue to assume, as in the previous section, that 
a z = ms, with m an arbitrary integer. The most general 
free energy consistent with these symmetries is 



K 



F = fd*r{-\(V-iAm 



1*1 



(55) 



where the coordinates have been rescaled to obtain an 
isotropic gradient term. The "vector potential" A repre- 
sents changes in the applied field 5H = SHf,y + H c z, with 
A x = 0, A y = qH c /H b , and A z = qSH b /H b . The form of 
this caunling follows from the transformation properties 

of <i>MeI 

Eq.|35] assumes a local form of the free energy. Ad- 
ditional non-local interactions arise due to interactions 
with long wavelength fluctuations in the density and tan- 
gent fields. The most relevant (near the critical point) of 
these couplings is 



9>— 5n 



"7 



d-W|$| 2 , 



(56) 



where the correlations of Sn are determined from Eqs.|2^ 
and 0. 

When SH = A = 7 = 0, Eq.||| is the free en- 
ergy of an XY model with an m-fold symmetry break- 
ing term. A second order freezing transition occurs 
within Landau theory when v > and r cx T — T s 
changes sign from positive (in the liquid) to negative 



8 



(in the smectic). The renormalization group (RG) scal- 
ing dimension, A m , of the symmetry breaking term is 
known experimentally in three, dimensions to be A m s» 
3 - 0.515m - 0.152ra(m - 1).E3 For m > m c w 3.41, the 
field <? is irrelevant (A m sC. 0), and the transition is in 
the XY universality class.ES The magnetie-.fields used by 
Kwok et. al.E3 correspond to m = 9 — 11, E-3 well into this 
regime. The static critical behavior is characterized by 
the correlation length exponent v « 0.671 ± 0.005 and 
algebraic decay of order parameter correlations at T s , 



(*(r)3*(0)> 



1 



,i+v ' 



(57) 



with 77 w 0.040 ± 0.003.E3 

To study the effects of coupling to long wavelength 
fluctuations when 7 ^ 0, we first satisfy Eq.Oby defining 
an auxiliary "displacement" -like field w via 



Sn = -V 



T = dyW. 



(58) 



After this change of variables, Eq.^?] becomes 

Fw = 2^1 / d3r | Cll l V ^ ' w ! 2 + c M,\\\d y w x \ 



+CiA,±\d y w z \ 2 



(59) 



where we have taken the q = limits of the elastic moduli 
to study the critical behavior, and dropped the Vjp term 
which only couples to w at finite q z . Eq.pffl then becomes 
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ePrV 1 • w|$| 



(60) 



Eq|30] is an anisotropic form of a coupling studied pre- 
viously in the context of the compressible Ising model, 
in which w describes the phonon modes of a compress- 
ible lattice on which the spins resideEj As shown in ap- 
pendix the techniques developed for that problem give 
the renormalization group eigenvalue A 7 = a/2i> for this 
coupling at the critical point. Since a = 2 — Zv —0.01 
is negative, the long wavelength density fluctuations are 
irrelevant for the critical behavior. 



VI. SMECTIC PHASE 

A. Static Behavior 

Deep in the ordered phase [r < 0), amplitude fluctu- 
ations of <E> are frozen out. Writing <i> = yJ\r\/ve 2 ' Klu / a , 
Eq.pq becomes, up to an additive constant, 



smectic 



J dh^(Vu-A) 2 - <7Cos27ra/sj 



(61) 



where a = ms, a = Air 2 \r\K / 'a 2 v, g — g(\r\/v) m / 2 , and 
the reduced vector potential is A = A/q. The displace- 
ment field u describes the deviations of the smectic layers 
from their uniform state. The sine-Gordon term is an 



effective periodic potential acting on these layers. Asia 
well known from the study of the roughening transition,tZI 
such a perturbation is always relevant in three dimen- 
sions. The smectic state is thus pinned at long distances 
(i.e. the displacements u of each smectic layer are local- 
ized in a single minima of the cosine). 

FIG. 5. Fluctuating density wave in the smectic state. Dis- 
placements of the layers from their mean positions are de- 
scribed by the field u. 

To further characterize the smectic phase, we consider 
the transverse magnetic susceptibility, which defines the 
macroscopic tilt modulus, 



-44, j 



(62) 



H c =0 



The field H c attempts to tilt the smectic layers. How- 
ever, A cx H c is an irrelevant operator in the smectic 
phase (as can be easily seen by replacing the periodic 
potential by a "mass" term a u 2 ). This implies that the 
smectic layers do not tilt under weak applied fields, i.e. 
d(d y u) / dH c \H c =Q = 0. Naively, this implies an infinite 
tilt modulus. 

A more careful treatment shows that 044.^ actually 
remains finite in the smectic phase. To compute 044.^ 
from first principles, we use the thermodynamic relation 



-Ait 



df 



(63) 



where / is the full free energy of the system, including 
a smooth part /o not involving $ and not included in 
Eqj55|, i.e. 



/ = /o - k B T In Z^>. 



(64) 



To evaluate Eq|33|, we need to consider in detail the de- 
pendence of the free energy and the coefficients in Eq.|55| 
on H c . This dependence arises in two ways, because 
an applied H c can be decomposed into a rotation and a 
scaling of the full field H. If the system were fully rota- 
tional invariant, the rotational part would enter F purely 
through the "gauge-invariant" coupling to A of Eq,55 
However, anisotropy breaks this invariance, leaving in 
stead only an inversion symmetry under z — > — z. The 
inversion symmetry-allows for a quadratic dependence of 
r and of f on H c r3 The scaling part also contributes 
quadratic dependence, which may be combined with the 
previous effect. Taking both into account, and matching 
the tilt modulus to the tilt modulus of the liquid phase 
(i.e. with <I> = 0) leads to 



S c = ^Im($*9 y $) 
rib 



+ C 



H, 



(65) 



where 044.^0 is the tilt modulus obtained from anisotropic 
GL theory (without accounting for the discreteness of the 
layers) and r" = d 2 r/dH 2 \H c =o- 

Eq.pq has a simple physical interpretation. The first 
term is the contribution to B c from tilting of the lay- 
ers (described by a phase shift of $). This term is 
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zero for small fields H c due to the cosine pinning po- 
tential. Even when the layers retain a fixed orienta- 
tion perpendicular to the c axis, however, the trans- 
verse field can penetrate via the second term. Such mo- 
tion arises microscopically from a non-zero equilibrium 
concentration of vortices with large kinks extending be- 
twee n nei ghboring smectic layers, as suggested in sec- 
tion IV . Eq.[35| predicts a non-divergent singularity 
c 44 ^(T) - c 44 ,x(T s ) ~ |T - Tsl 1 -" at the critical point, 
where a is the specific heat exponent. 

At low temperatures in the smectic phase, we can esti- 
mate the tilt modulus in terms of properties of kinks. In 
zero field, the concentrations of large kinks carrying mag- 
netic field in the +z and z directions are equal, leading 
to zero net field along the c direction. For H c ^ 0, the 
energy of a kink depends upon its orientation due to the 
—B c H c /Att term in the GL free energy, yielding 



E± w £ik ± ms(f>oH c /ATT. 



(66) 



The difference in the concentrations of up and down kinks 
takes the activated form 



71+ — Tl- 



B & /,, rp I ms<boH r 
e -*WfcBT sinh VO c 



(67) 



where wik ~ y^ej_/U p ms, estimated from Eq|25]with , 
ms. Since B c = (n + — n )ms(j)o, Eq.p2| yields 



C44,_L 



y/e ± /U p k B T 



t>oms 



,E lk /k B T 



(68) 



i.e. a large but finite tilt modulus. 



B. Dynamical Behavior 

Very similar phenomena occur in the dynamics of the 
smectic phase. To study them, we need the equation of 
motion for $. On the basis of symmetry and the lack 
of obvious conservation law*, a natural conjecture is that 
of overdamped "model A"tJ dynamics. Indeed, a care- 
ful treatment using the general formalism of section II B 
gives (see appendix |5|) 



■SF a 



5<f>* 



ip,J x $ - ij, 



(69) 



where p — q(ftono/c and r?(k) = inoqrj z (qz + kl Eq.|69] 
is remarkably similar to the model E dynamics^] for the 
complex "superfluid" order parameter <E>, where now J x 
plays the— role of the "electric field" in the Josephson 
couplingM3 The actual electric field is £ x = j v , z 4>o/ c : 
leading via Eq.[l5] to (see appendix 



"000 

2qc 



Im($*9 t $) 



B 



+ (1 - |$| 2 /2) ( — ) p xx , n J x , 



(70) 



where p xx ^ n is the normal state resistivity in the x direc- 
tion, whose appearance in the last term follows from the 
relation (n </>o/c) 2 /7BS ~ (B / H c2 )p xx ,„. 



FIG. 6. Sliding of a kink (thick curved line) along the field 
direction, viewed along the x axis for the case m — 2. Dashed 
lines indicate the copper-oxide layers. As the kink moves 
along the y axis, net vorticity is transported in the z direction. 
Such motion produces a finite resistivity in the smectic phase. 

Eq.[70| is interpreted in close analogy with Eq.^35|. In 
the absence of pinning due to the periodic potential in 
Eq.^l|, an applied force induces a uniform translation of 
the layers, and thus a net transport of vortices. In the 
order ed p hase, where $ = y / \r\/ve 27Tm ^ a , the first term 
in Eq.[70| becomes proportional to the velocity dtu. The 
second term contributes even when the layers are con- 
stant. It results microscopically from the motion of equi- 
librium vortex kinks, which can slide unimpeded along 
the y axis and thereby transport vorticity along the z axis 
(see Fig.^J) . Such flow at "constant structure" is analp-. 
gous to the permeation mode in smectic liquid crystals." 

The presence of this defective motion implies a small 
but non-zero resistivity at the L-S transition. Near T s , 
Eq.|7^ predicts a singular decrease of the form p xx (T) — 
p xx (T s ) ~ \T S — T| 1_Q , similar to the behavior of the tilt 
modulus. At lower temperatures (but still within the S 
phase) transport occurs via two channels. The perme- 
ation mode gives an exponentially small linear resistivity 
p xx ~ exp(— Eik/ksT) (above T s , single layer kinks give 
p xx - exp(-£ k /fc B T), with E k w E^/m). 

Non-linear transport occurs in parallel to the above 
linear processes, via thermally activated liberation of vor- 
tex droplets, inside which u (or u z in the crystal phase) 
is shifted by s (see FigJ^). Such a droplet costs a surface 
energy, due to the creation of a domain wall between 
smectic regions shifted by s. The domain wall surface 
tension <7o is estimated from 



Co 



K / S V 

2 \w) 



gw, 



(71) 



where w is the width of the domain wall. The first term 
represents the elastic cost of the shift in u, while the 
second is the pinning energy. Minimizing Eq.|7l] gives 
w ~ y/ K /s s an d Co ~ y/ngs. This surface energy must 
be balanced against the Lorentz force in the interior, so 
that the energy of a droplet of linear size L is 



-^droplet ~ CqL 2 



JBs 



(72) 



Eq.[72| gives a critical droplet size L c ~ ifivyc/ \JB) 
and an energy barrier Eb ~ (c/JB) 2 (kj) 3 ^ 2 s. Thermal 
activation therefore gives 



(73) 



where J c - (c/ B)(Kg) 3 / 4 (s/k B T) 1 / 2 . Similar non- 
linear IV relationahave been obtained previously for vor- 
tex/Bose glasses,aQ but our result is more closely related 



to surface mobi. 
crystal surfaces.E 



y below the roughening transition on 
] Unlike these proposed glass phases, 
the smectic should always exhibit a nonzero linear resis- 
tivity as J — > 0. 
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FIG. 7. A two dimensional cut through a droplet configu- 
ration of the smectic layers, drawn for the case m = 2. The 
full three-dimensional droplet has a spherical droplet. Inside 
the droplet (outlined in gray), the layers are shifted by one 
CuC>2 double layer spacing, u — > u + s. 



VII. SUPERSOLID ORDER AND THE SMECTIC 
TO CRYSTAL TRANSITION 

A. Supersolid Nature of the Smectic Phase 



In sections VIA and VI B, we have seen that the re- 



sponse functions in the smectic phase retain many of 
the features of the vortex liquid. Both the tilt modulus 
and conductivity remain finite, despite the pinning of the 
smectic density wave by the CuC-2 layers. As discussed 
earlier, both phenomena are explained by the existence 
of an equilibrium concentration of large vortex kinks ex- 
tending between successive occupied vortex layers. These 
kinks facilitate both transverse magnetic penetration and 
dissipation for currents along the x axis. 

This behavior is strikingly similar to the picture of "su- 
persolid" vortex arrays recently proposed in Ref. for 
fields parallel to the c-axis. In the supersolid, a finite 
concentration of interstitials or vortices are present in the 
vortex lattice, and both the tilt modulus and conductiv- 
ity in the presence of weak pinning are finite. Such a 
supersolid phase is distinct from the Abrikosov solid in 
that it supportslong range crystalline order coexists with 
a finite expectation value of the boson order parameter 
ip, i.e. 



(V>(r)V>*(0)) -» Const. 



(74) 



as r — > oo. The supersolid must occur at sufficiently high 
magnetic fields, but its [existence elsewhere in the phase 
diagram seems unlikelyJia 

Using this characterization of broken U(l) symmetry 
(under ip — > ipe 1 ®), the vortex smectic is always in a su- 
persolid phase. As in Ref. |l^, this can be seen by consid- 
ering the correlation function of tp's. Note that ip(r) de- 
stroys a vortex line at position r are "0*( r ) creates a line in 
the coherent state path integral formalism. Because there 
is always a finite probability of finding a kink connecting 
the points and r, Eq.74 is indeed satisfied. With this 
understanding, the second terms in Eqs.|7^ and ^35] have 
an additional complementary interpretation. They corre- 
spond to the contributions from the "superfluid fraction" 
of a two-fluid system with "superfluid" (kink) and "nor- 
mal" (smectic) parts. 

In addition, the concept of symmetry breaking im- 
plies that a continuous transition from the flux liquid 
state (with (ip) ^ 0) to a smectic (translationally or- 
dered in one direction) phase must necessarily retain su- 
persolid order. For the smectic phase to appear with 
(ijj) = would require simultaneous breaking of the dis- 
crete translation group and restoration of the E/(l) sym- 
metry. Such a double critical point can only occur by 
tuning two parameters (one in addition to the temper- 
ature) or through a first order transition. The physical 



arguments of sections IV C and VIA-VIB, of course, im- 
ply the stronger condition that the smectic phase must 
be supersolid at all temperatures. 



B. Consequences for Further Transitions at Low 
Temperatures 

At lower temperatures, provided point disorder re- 
mains negligible, the vortices will order along the x axis 
as wellMJ What is the nature of this two-dimensionally 
ordered phase? 

The different possibilities may be classified by the or- 
der in which the symmetries are broken. At the lowest 
temperatures, we expect the system to prefer a true solid 
phase, with broken translational order in both directions 
(in particular (pii) 7^ 0, where pii is the amplitude for 
periodic density variations along the a;-axis. Recall that 
p± = $ is the amplitude for density waves along, z.) and 
a restored U (1) symmetry (i.e. no interstitials) J13 To con- 
nect this state with the smectic phase in which (p\\) = 
and (ip) 7^ requires two changes of symmetry. 

Three monotonic choices of symmetry breaking are 
shown in Fig.|| In the scenario (a), upon lowering the 
temperature from the smectic phase first the U(l) sym- 
metry is restored, and the translational symmetry along 
the x axis is broken at a lower temperature. As remarked 
in the previous section, however, the intermediate non- 
supersolid smectic phase that appears in this sequence is 
impossible, so this sequence cannot occur. 

Two physical choices remain. The smectic may go di- 
rectly to the normal solid in a first order transition which 
breaks the translational symmetry and restores the U{1) 
invariance simultaneously, as shown in Figj8|(b). The last 
possibility, illustrated in Fig.||(c), is that of an inter- 
mediate supersolid phase between the smectic and the 
interstitial-free solid. In this case both the low tem- 
perature phase transitions may be second order. The 
supersolid-solid critical behavior is described in Ref. [ll]. 

FIG. 8. Three routes of symmetry breaking at low tem- 
peratures connecting the smectic phase to the vortex solid. 
Choice (a) is ruled out on physical grounds. Note that 
(p±) = (<&) 7^ in all cases. 

The smectic-supersolid transition is once again a freez- 
ing transition at a single wavevector, and is potentially 
describable by a Landau theory like Eqj5^. Because the 
modulating effect of the underlying crystal lattice is much 
weaker in the x direction, we expect g sa is a good ap- 
proximation in this case, which leads to pure XY behav- 



VIII. INCOMMENSURATE PHASES 

As is-.well known from the study of the sine-Gordon 
model,ll3 a large incommensurability can be compensated 
for by energetically favorable "solitons" , or walls across 
which u — > u + s (see Fig j^) . Solitons begin to proliferate 
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when their field energy per unit area (Xfi c id ~ —kAs ex- 
ceeds their cost at zero field, ctq ~ \/k9 s (estimated from 

Physically, these solitons correspond to extra/missing 
flux line layers and walls of aligned "jogs" for 5H along 
the b and c axes, respectively (see Figj9j) . In the for- 
mer case, this leads to an incommensurate smectic (IS) 
phase, whose periodicity is no longer a simple multiple 
of s. For <5H || z, the solitons induce an additional peri- 
odicity along the y axis. This tilted smectic (TS) phase 
has long range translational order in two directions .EJ 
The analogous tilted crystal (TX) phase is qualitatively 
similar, but has long range order in 3 directions. 

For larger H c , as the angle between the field and 
the CuC>2 layers becomes large, intrinsic pinning and 
anisotropy no longer favor the smectic state. As shown in 
Fig.^j, we therefore expect the L-TS and TS-TX phase 
boundaries to merge in this regime. The direct L-TX 
transition is necessarily first order. 

FIG. 9. Kinked configurations (neglecting fluctuations and 
drawn viewed along the x axis) of the smectic layers for mag- 
netic field perturbations (a) along the b axis, and (b) along 
the c axis. 

In conventional CITs, entropic contributions generate 
additional interactions between domain walls which ac- 
tually dominate over the bare energetic repulsions when 
the inter-soliton spacing I — > oo. To estimate their mag- 
nitude here, we use the well known logarithmic roughness 
of a 2d interface 

((Mx)-M0)) 2 )~mM, (75) 

where h is the height of the interface and the coordinate x 
parameterizes its position in the base plane. For solitons 
spaced by £, collisions between neighbors generally occur 
only once h > £, so that the size of roughly independently 
fluctuating regions x ~ exp(^ 2 ). The entropy loss due 
to this constraint scales with the number of collisions 
(L/x) 2 , so that the areal free energy cost per wall is 

/ col l. = -TAs coll . ~ Te~ e . (76) 

Since the energetic interactions in the smectic scale expo- 
nentially (like exp(— £/w)) at long distances, the collision 
free energy is actually negligible as £~=r* oo, unlike the 
situation for lines in 1 + 1 dimensions £3 The free energy 
density in the incommensurate phases is thus 

f I "' _i_ ^ c -i/w / 77 \ 
/soiiton ~ j- + —e 1 , (77) 

where a = crfidd + <Jq < is the total areal free energy of 
the soiiton; A and w set the energy and length scales of 
the soiiton interactions. At low temperatures we expect 
w ~ A and A ~ £o/a x , while near T s , Eqj6l] gives w ~ 
y/ n/gs (c.f. Eq.fnj) and A ~ gl. Minimizing Eq.|77] gives 
a soiiton separation I ~ wln(A/|cr|) near the CIT. 

In the TS phase, net vortex motion along the c axis 
occurs by sliding soiiton walls along the b direction. 
The resulting electric field is proportional to J and the 



soiiton density, leading to an additional contribution to 

the resistivity which vanishes at the CIT like p x ° x ton ~ 
^ oliton /ln(A/|cr|). 

A single soiiton wall in the IS phase, because it is par- 
allel to the Cu02 layers, experiences a periodic potential 
along z. From studies of the roughening transitionj£3 it 
is known that such a periodically pinned wall may be in 
either a rough or smooth phase. If the walls are individu- 
ally smooth, thermal fluctuations are negligible, and the 
assembly of solitons is well described by an effectively, 
one-dimensional elastic chain in a periodic potential.Ej 
Because they are pinned separately into minima of the 
potential, they do not contribute to p xx . If they are 
rough, they wander logarithmically and eventually inter- 
act with their neighbors. The appropriate coarse-grained 
description beyond this interacting length scale is an elas- 
tic stack of domain walls. The configuration of such a 
stack is described by a second displacement field «d w , 
with a free energy of the same form as Eq.[n] (but with 
different values of k and g). The statistical mechanics for 
the ttdw field is thus equivalent to that of the original u 
variable. The preceding analysis must then be repeated 
within the new effective free energy. 

Because of the aforementioned complexity of the one- 
dimensional problem, we have not attempted to deter- 
mine the true long distance behavior of the soiiton array 
in the IS phase. Because the permeation mode in the 
commensurate smectic already provides a finite tilt mod- 
ulus and non-zero resistivity, however, we expect that 
these more subtle effects will have only weak experimen- 
tal implications. 

As the temperature is increased within the IS or TS 
phases, the system melts into the liquid. To study such 
transitions, we perform the dilation $ — > $exp(iA- r). 
Only the g term is not invariant under such a gauge- 
like transformation. It becomes oscillatory and therefore 
does not contribute to the critical behavior at long wave- 
lengths. The IS-L and TS-L phase transitions are thus 
XY-like. 

The shape of the CIT phase boundary is of particular 
experimental interest. In the mean field regime, this is 
obtained from the condition a = as 5H ~ |r| T , with 
Tmf = (m— 2)/4. By the usual Ginzburg criterion, mean 
field theory breaks down for \r\ < (k B Tv/ K 3 ^ 2 ) 2 . To 
determine the shape of the phase boundary in this critical 
regime, we follow the RG flows out of the critical region 
and repeat the preceding analysis with the renormalizcd 
couplings determined by matching when |r| is order one. 
Then SH R - £ Xh 5H and g R ~ £ Am #, with £ - \r\-' y . 
Rotational invariance at the rescaled fixed point (g = 0) 
implies that the field exponent is exactly Ah = 1 (see 
appendix ^). Using these renormalized quantities, we 
find 

Tcrft. = (|A m |+ 2W2«4.9 -7.2, (78) 

for the fields used in Ref. [l^. 

The IS-L and TS-L phase boundaries are non-singular 
and are determined locally by the smooth 5H dependence 
of r. In particular, for small H Cl the TS-L phase critical 
temperature is 
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T(H C ) =T„- —Hi 



(79) 



where r' = dr/dT\ 



t=t 3 ,sh=o- 



IX. INFLUENCE OF DISORDER 

Lastly, we consider the effects of weak point disorder, 
which couples to the density of vortices according to 



F d = J d\V d {v)n{v), 



(80) 



where V d (r) is a random potential, which, for point im 
purities, is short range correlated in space and narrowly 
(e.g. Gaussian) distributed at each point. Using Eq.54 
F d can be rewritten, up to less relevant terms, as 



F d = J dW d (r)Re{$(r)e i<? *}, 



(81) 



in the smectic phase and in the liquid sufficiently near 
T s . To bring Eqj8l] into a more standard form, we define 
a complex random field V d = V d e lqz , in terms of which 



(82) 



Because of the oscillatory e lqz factor, V d and Vj are es- 
sentially uncorrelated at long wavelengths. Eq ]82] i s the 
simplest "random field" XY perturbation of Eq!p5L and 
the resulting model is known in statistical mechanics as a 
random field XY model with an m-fold symmetry break- 
ing term. 

Before discussing the critical behavior of such a theory, 
it is natural to consider the effect upo n the o rdered state. 
In the smectic phase, using $ = y/\r\/v exp(27rra/a), 
Eq.^2| becomes 

F d = J d\^\j^ (y d ¥ e 2mla + V d e~ 27Ilu/a ^j . (83) 

If the disorder is weak relative to the periodic potential 
(i.e. \V d \ <C g), it is naively justified to replace the cosine 
in Eq.|6l| by the "mass" term (27r / ' s) 2 gu 2 / '2. Such a mass 
term gives a large penalty for excursions of the layers 
with u s, so the randomness in F d appears irrelevant. 

By ignoring the periodicity of the cosine, the above 
approach does not consider the possibility of disorder- 
induced solitons. To study the stability of the smectic 
to such topological defects, consider a region of size L in 
which the displacement field u is shifted by s, so as not 
to incur any bulk energy cost from the intrinsic pinning. 
Within this region, F d contributes an energy of random 
sign of order \V d \L d / 2 in d dimensions. On the boundary 
of the region, however, the cosine does contribute, costing 
an energy ~ gU 1 ^ 1 . For d > 2, the boundary energy 
grows more rapidly with L, and the net energy is always 
positive, provided g > |Vd|s 1_d / 2 . Thus we see that the 
smectic phase remains stable to weak disorder, even once 
solitons are taken into account. 



Note that this result is in strong contrast to the 
Larkin-Ovchinikov argument that the .Abrikosov lattice 
is unstable to arbitrarily weak pinning.ll3 Physically, the 
instability is prevented, at least for weak disorder, by the 
periodic pinning potential which increases the stiffness of 
the smectic displacement field. The result can be under- 
stood, however, on more general symmetry grounds. For 
non-zero g, the system does not have a true continuous 
translational symmetry in the z direction, but only the 
discrete symmetry under translations by s. Because the 
symmetry is discrete, there is no Goldstone mode in the 
ordered (smectic phase) - i.e. phonons are massive. It is 
now well known that for random field models with dis- 
crete symmetries (e.g. the random field Ising model, to 
which our model corresponds when ttVt ^ ne ordered 
phase survives above two dimensions E3 Indeed, the ar- 
gument given above for stability against droplet solitons 
is a restatement of the Imry-Ma argument first used for 
the random field Ising modelo 

In the incommensurate (IS and TS) phases, where the 
periodic pinning g is effectively zero, the Imry-Ma ar- 
gument no longer applies. In these phases, the original 
Larkin-Ovchinikov picture holds, and the distortions in 
|u(r)| 2 must grow on long length scales, destroying the 
long range translational order of the layers. The nature 
of the resulting phase is unclear: it may be a "smectic 
glass" , analogous to the proposed vortex glass phase for 
more isotropic systems, or it may simply be a strongly 
correlated liquid, with slow relaxation times. The same 
considerations hold for the more ordered phases at low 
temperatures, since the translational order along the x 
axis lacks the intrinsic pinning required to prevent the 
Larkin-Ovchinikov instability. 

Turning to the critical behavior of the L-S transition, 
the analysis becomes more subtle. The random field per- 
turbation in Eqj8| is a relevant perturbation at the XY 
fixed point (and indeed at any 0(n) fixed point), so the 
critical behavior is certainly altered. Naively, a 6 — e 
expansion may be made for the critical behavior of the 
0(n) randomiield model, which has a zero temperature 
critical point o Within such a perturbative expansion 
near 6 dimensions, the symmetry-breaking term appears 
irrelevant. There are two potential problems with this 
approach. Firstly, if the symmetry-breaking term indeed 
remains irrelevant at the new critical point, it would be 
an example of a three dimensional random field XY criti- 
cal point. The random field XY model, however, because 
it has a continuous symmetry, does not even have a stable 
ordered phase in three dimensions. Although there may 
not be an obvious contradiction involved in this scenario, 
the physical meaning is certainly unclear. One possible 
resolution is that the symmetry-breaking term becomes 
relevant at some higher dimension (greater than 4) . The 
second problem is the status of the 6 — e expansion itself, 
which has been proven to break down, at least via non- 
perturbative corrections (but possibly moije-^trongly) for 
the case of the random field Ising model.cll The consis- 
tency of the 6 — e expansion, even perturbatively, has not 
yet been determined. Regardless of the success or fail- 
ure of this theoretical approach, experimental work on 
random field Ising systems has demonstrated the subtle 
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types of behavior possible for such zero temperatur e crit - 
ical points. Fortunately, as discussed in sections VI A - 
VII , the supersolid nature of the smectic (with the rela- 
tively fast permeation mode) implies that slow dynamics 
for the smectic ordering will not have a strong impact on 
transport and magnetization experiments. 

Finally, consider the L-IS and L-TS transitions in the 
presence of disorder. Disorder strongly effects the be- 
havior at such CITs, because it modifies the wandering 
of a single domain walloo Fortunately, the effects of 
randorrv-4ji&l|d disorder on a single interface are known 
exactlySH In contrast to Eq.[75], the height fluctuations 
of the interface grow like 



(ft(x) - h(o)y 



l 2 C 



(84) 



where C = (5 — d)/3 — 2/3 in three dimensions. Also un- 
like the pure interface, the free energy fluctuations within 
a region grow with length scale, so that the cost per col- 
lision scales like \x\ e , where 6> = <i-3 + 2C = 4/3. The 
areal collision free energy per wall is thus 



coll. 



2-0) /( 



i/e. 



(85) 



This dominates over exponential energetic interactions 
for large I. The full soliton free energy is thus 



soliton 



i 



Ad 

(2 ' 



(86) 



where measures the strength of the disorder-induced 
collision interactions. Minimization of Eq||^ gives t ~ 

l/kl. 



X. CONCLUSIONS AND APPLICATIONS TO 
HELIUM FILMS 

We have studied the behavior of vortex arrays sub- 
jected to a one-dimensional periodic potential transverse 
to the magnetic field. Such a potential, which is induced 
by the layered structure of the high temperature cop- 
per oxide superconductors for fields oriented in the a-b 
plane, favors an intermediate smectic phase between the 
vortex lattice and flux liquid. The commensurate smec- 
tic state is supersolid, and has in consequence a nonzero 
finite resistivity and tilt modulus, despite being pinned 
by the periodic potential. Including incommensurability 
effects leads to the rich phase diagrams of FigJ|. The 
experimental signature of the smectic is the appearance 
of Bragg peaks in the structure function along a single 
ordering axis interleaving the trivial peaks induced by 
the layering. We also expect a greatly reduced resistivity 
for currents transverse to both the layering and magnetic 
field axes, and a cusped phase boundary describing the 
response to small fields perpendicular to the layering di- 
rection. The qualitative behavior of the tilt modulus and 
resistivity near the liquid-to-smectic transition in a field 
aligned perfectly with the ab-plane is shown in Fig. 10 



FIG. 10. Tilt modulus and resistivity at T s for a perfectly 
aligned commensurate field. Not to scale: the resistivity could 
appear to drop to zero and the tilt modulus could appear to 
diverge to infinity in all but the most precise experiments. 
The open circles denote \T — T s | 1_a singularities, where a is 
the specific heat exponent. 



Using the boson mapping,H€3 these results can be 
extended to real two dimensional quantum mechanical 
bosons at zero temperature. The analogous quantum 
smectic phase might be studied in helium on a periodi- 
cally ruled substrate. Such a substrate might be approx- 
imated by crystalline facets exposing a periodic array of 
rectangular unit cells with large aspect ratio. There are, 
however, a number of difficulties inherent in this exten- 
sion. In particular, the interaction between helium atoms 
is not purely repulsive; it is reasonably well described by 
a Lennard-Jones potential with a minimum at an inter- 
atomic spacing of a few Angstroms. To obtain a sub- 
strate with a small enough period to affect the physics 
on these length scales is an experimental challenge. Were 
these interactions purely repulsive, one could probably 
overcome this difficulty by working with a dilute system. 
With an attractive tail to the potential, however, a low 
density helium film would likely phase separate into he- 
lium rich and helium poor regions, making intermediate 
densities inaccessible. Appropriate experimental condi- 
tions may nevertheless be achievable for small values of 
m (the number of periods of the potential per period of 
the smectic density wave). The depth of the minimum in 
the effective pair potential, moreover, could be reduced 
somewhat by a careful choice of substrate. 

The case m = 2 has been explored numerically in a 
Bose Hubbard model in Ref. |5^. These authors indeed 
find a smectic phase, which they denote a "striped solid" , 
with order in reciprocal space at q = (tt, 0). Their results 
for the structure function and superfluid density appear 
to be in good agreement with the predictions of section 
[Vl| (see in particular their Fig.ll. The superfluid density 
of the boson system maps onto the inverse tilt modulus 
c^l of the flux lines .0). Unfortunately, a detailed com- 
parison of the critical behavior with the theory is beyond 
the resolution of the available data. 

At a more general level, the liquid-smectic transition 
treated here appears to be the only known case of con- 
tinuous quantum freezing in 2 + 1 dimensions. The smec- 
tic phase, moreover, is perhaps the simplest example of 
a quantum phase intermediate between solid and liq- 
uid. The techniques developed here may be useful in un- 
derstanding other quantum phases of mixed liquid/solid 
character. One particularly intriguing example is the 
"Hall solid" proposed in Ref. |33|. Through the Chern- 
Simons mapping,L3 it can be related to a supersolkLphase 
of composite bosons (electrons plus flux tubes) £3 This 
phase, an analogous "Hall smectic" and "Hall hexatic," 
and the modifications of the current theory to account for 
the long-range Coulomb and Chern-Simons interactions 
are discussed in Ref. B5l 
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APPENDIX A: DERIVATION OF 
HYDRODYNAMIC EQUATIONS 

The continuity equation (Eqjl3|) follows from differen- 
tiation of EqsJ|. For example, 

h = — [S'(x — Xi)S(z — Zi)ii + 5(x — Xi)S'(z — Zi)zi\ 

i 

= -V ± • Yj 6 ^ 1 - ~ r -Li) f "-Li = -V_L • j„, (Al) 

i 

where 

J" = £ S ( T± ~ r -i-t) f '-L*- ( A2 ) 

i 

Eq.|l4| is derived analogously, giving the tangent current 
tensor 



3p 



a = y^S(r± - r_u) 



dxf dxf dxf dxf 
dt dy dt dy 



(A3) 



Eqs. A2 and A3 are completely general, and do not de- 



pend upon the detailed dynamics of the vortex system. 
This additional physics is included in the constitutive 
equations (EqsJlJ and |l6|). To derive them, we need the 
equation of motion, EqW Inserting this into Eq. A2 gives 



Pk 



SF 



Sr±i(y) 



nf . 



(A4) 



The functional derivative with respect to r±i can be 
transformed via the chain rule 



SF = / d J r 



" SF . , , SF . , ' 
Mr) + Srjr) " 5T{T) . 
where the variations 5n and St are 

Sn(r) = -V_l • <5(rj_ - r±i)Sr± t (y) 



(A5) 



8t(t) = -d a ^S(r± - r_u) 

i 



dr± 

dy 

d5rj_j 
dy 



■Sxf(y) 



(A6) 



Substituting Eqs.|Aj and |A6| into Eq|A4j then gives Eq|15 
The constitutive equation for jp a is derived analogously, 
by premultiplying EqjT] with dxf /dy and carrying out 
the same steps as before. 



APPENDIX B: LANDAU THEORY OF THE 
SMECTIC TO CRYSTAL TRANSITION 

We assume that H is in the ab-plane and commen- 
surate smectic order is already well established, and ask 
how a two-dimensional vortex modulation then arises at 
a lower temperature. The modulated vortex density now 
takes the form 



n(r) = n Re< 1 + $(r)e~ iqz + V>i(r)e" 



-»GlT 



-iG 2 r 



(Bl) 



where $(r) is the (large) smectic order parameter, and 
Gi and G2 are reciprocal lattice vectors lying in the (x — 
z) plane with G\ x = —G2 X 7^ satisfying 

qz + Gi + G 2 = 0. (B2) 

The six vectors ±gz, ±Gi, and ±G 2 form a distorted 
hexagon of minimal reciprocal lattice vectors. All other 
reciprocal lattice vectors in the crystalline phase are lin- 
ear combinations of this set, which reflects an anisotropic 
vortex lattice in real space. The corresponding set of re- 
ciprocal lattice vectors for a square lattice is illustrated 
in Fig§ 

The complex amplitudes ipi(r) and ip2 (r) are small 
near the transition, and the Landau free energy differ- 
ence SF between the smectic and crystalline phases takes 
the form 



SF 



d*r 



y |VV^i| 2 + y |VV> 2 | 2 +gV^ ■ VV2 



+ o(l^l| 2 + l^l 2 ) + w($VlV>2 + + 



(B3) 



We have equated the coefficients of gradients in all three 
directions for simplicity. Within mean field theory, crys- 
talline order can arise via a continuous phase transition 
whenever r < 0. The neglected higher order terms fix 
the magnitudes of ipi and tp2, I "01 1 = IV^I = i>o below the 
mean-field transition temperature. To study the true 
transition (which occurs for r = r c < due to thermal 
fluctuations), we set 

^ = ^ e iei(r) , V2=V'oe ie2(r) . (B4) 

Upon neglecting a constant, the free energy becomes 



SF 



K, 



IV0J 



K 



\V9 2 



+gV9x ■ V6 2 +wcos(9 1 - 6 2 ) 



(B5) 



where K — Ki/jq, g — gipQ, w — w\^\i^, and we have 
assumed the phase of the smectic order parameter $ is 
locked to zero by the periodic pinning potential. Eq.B5 
represents two coupled XY models with phases locked 
by the cosine. This term forces 9± w 0% = 0, and the 
phase transition falls in the universality class of a three- 
dimensional XY model with effective free energy 



SF X Y ^(K + g) d 3 r\V9 



(B6) 
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APPENDIX C: EFFECT OF LONG 
WAVELENGTH FLUCTUATIONS AT THE S-L 
CRITICAL POINT 

To determine A 7 , we assume the usual scaling form of 
the free energy at the critical point, 



/(r,7)=rt(7^), 

where g is an unkn own scaling function. Using £ 
and 2 — a = dv, Eq.Cl gives 



9V 
$7 2 



g(0) ~r 



2 — Oi— 2\-y 



(CI) 



(C2) 



7=0 



The same quantity can be calculated directly, however, 
by differentiating the partition function to obtain 



$7 2 



= [ d 3 x(V ± -w(x)V ± -w(0))(S(x)^(0)), 
=0 J 

(C3) 



where the "energy" operator E{x) = |$(x)| 2 . The angu- 
lar brackets indicate expectation values in the decoupled 
theories. The energy-energy correlations take the scaling 
form 



(E(x)E(0))~r^-^h(\x\/0, 



(C4) 



where h(x) i s another unknown scaling function. The 
long-wavelength fluctuations, determined from Eq.59 
are 

(V_l-w(x)V±-w(0)) = nlk B T x 



C44,±K + c 44,||<7z 



e lqx . (C5) 

(27r) 3 cii(c 44 ^g 2 + c M ,\\Qz) + C44,±c 4 4,||9y 



Inserting Eqs.C4 and C5 into Eq.C3 and changing vari- 
ables x — > £x gives the scaling 



d"f 2 



,2(1- 



(C6) 



7=0 



Comparison with Eq.C2 then gives the desired result 
A 7 = a/2i/. 



APPENDIX D: EQUATIONS OF CRITICAL 
DYNAMICS 



To derive the critical equations of motion, we assume 
that important fluctuations occur only near q = and 
q = qz. Including the q = modes in Eqj54 gives 



n(r) = n [l + Re$e~ ig2 )] + Sn. 



(Dl) 



It is sufficient to keep only the small q parts of the tan- 
gent fieldES The finite q z modulation of the density in- 
duces a modulation of the vortex current, 



Jr 



(D2) 



where j VjU and j„ )S are the (slowly varying) uniform and 
smectic (i.e. per iodic) components of the vortex current. 
Inserting Ecj jDl| into EqJlJ and isolating the parts pro- 
portional to e~ lqz 



-4 1 



gives 
Sn 



+n <S>d y — 



no 
SF 

^7 



(v 



iqz) 



SF 



n M + nofthormaie 192 



no'I'Vj 



SF 
S(Sn) 

(D3) 



where the last term is included because the white noise 
fthermai has Fourier components at all wavevectors. Sim- 
ilarly, from the uniform component, one finds 



rj^,„ = -(n + Sn) 



-2Re 



#(Vj 



SF 

S{Sn) 



St 



_ SF 

~T a V ± — 

Sr a 



+ (n + Sn)f. (D4) 



Using Eqs.D2, the divergence of the current is 

V_l • j v = Vj. • Jv,u + Rc(V ± - iqz,)j v , s e~ iqz , 
which leads to the two continuity equations 



n d t $ - 



d t Sn 
(V, 



iqz) -jv,s 



o. 
0. 



(D5) 



(D6) 
(D7) 



Eqs.D2-D7 completely specify the dynamic s of the 
de nsity fluctuations Sn and <I>. Inserting Eq.D3 into 
Eq.D7, and neglecting irrelevant couplings gives the 



"model E" -like model of Eq. 69 



Eq.D4 can, in principle, be used to study the effects 



of the smectic ordering on the long-wavelength modes. 
The most physical application, however, is to determine 
the contribution of the smectic degrees of freedom to 
the electric field. The average bulk (i.e. q = 0) field 
is (£ x ) — (j v ,u,z)4>ol c - I n the s piri t of the Landau expan- 
sion, q = component of Eq. D4 can be reasonably ap- 



proximated by dropping terms with gradients of $ (the 
leading contributions from the first term involving Sn 
and t vanished automatically at q = = even without this 
approximation) , yielding 



r(£ c ) w n 4>ofz/c + 



iqcf)o 



SF 
5§* 



-SF 
$ — 



(D8) 



Recognizing the var iation of the free energy on the right 
hand side of Eq.D8, the substitution (c.f. Eq.|69|) 



SF 
S$* 



7b s 



ft*- 



4g 2 



*, 



(D9) 



leads immediately to Eq.[70|. 



APPENDIX E: EIGENVALUE OF A 

The eigenvalue Ah is determined completely by rota- 
tional invariance. To see this, consider Eq|5^ at the fixed 
point, i.e. with g = - 0. By making the transformation 
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<J> — * <£>cxp(iA • r), the term involving A may be com- 
pletely removed from F. The free energy is thus inde- 
pendent of A. After integrating out short-wavelength 
modes, no dependence on A can appear, so precisely 
the same transformation must eliminate the field depen- 
dence in the renormalized free energy. This operation is 
$ — » <f>exp(iA ■ r£), using the rescaled r, so the renor- 
malized vector potential must be 



Ar = £A, 



(El) 



which implies Ah = 1. 
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